Evidence for Correlations Between Nuclear Decay Rates and Earth-Sun Distance 
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Unexplained periodic fluctuations in the decay rates of '^^Si and '^^^Ra have been reported 
by groups at Brookhaven National Laboratory (^^Si), and at the Physikalisch-Technische- 
Bundesandstalt in Germany (^'^^Ra). We show from an analysis of the raw data in these experiments 
that the observed fluctuations are strongly correlated in time, not only with each other, but also with 
the distance between the Earth and the Sun. Some implications of these results are also discussed, 
including the suggestion that discrepancies in published half-life determinations for these and other 
nuclides may be attributable in part to differences in solar activity during the course of the various 
experiments, or to seasonal variations in fundamental constants. 
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FoUowing the discovery of radioactivity by Becquerel 
in 1896 [l|] an intense effort was mounted to ascer- 
tain whether the decay rates of nuclides could be af- 
fected by external influences including temperature, pres- 
sure, chemical composition, concentration, and magnetic 
fields. By 1930, Rutherford, Chadwick, and Ellis [3, 
p. 167] concluded that "The rate of transformation of 
an element has been found to be a constant under all 
conditions." (For decays resulting from K-capture, or for 
beta-decays in strong ambient electromagnetic fields, the 
situation is slightly more complicated, since these decays 
are infiuenced by the electron wave functions which can 
be affected by external pressure or fields d, 0, For 
•^^Si and ^^^Ra, which decay by beta- and alpha-emission, 
respectively, fluctuations in the counting rates (in the 
absence of strong external electromagnetic fields) should 
thus be uncorrelated with any external time-dependent 
signal, as well as with each other. In what follows we 
show that neither of these expectations is realized in data 
we have analyzed for ^^Si and ^^^Ra, thus suggesting that 
these decays are in fact being modulated by an external 
influence. 

Between 1982 and 1986, Alburger, et al. Q mea- 
sured the half- life of ^^Si at Brookhaven National Labo- 
ratory (BNL) via a direct measurement of the counting 
rate as a function of time. If N{t) denotes the num- 
ber of surviving atoms starting from an initial popula- 
tion A'o at i = 0, then the familiar exponential decay 
law, N(t) = Noe-^\ leads to iV = dN/dt = -AA^oe"^* 



Nit) 



as a function 



where A = ln(2)/ri/2. A plot of In 
of time is then a straight line whose slope is A, which 
then gives the half- life Ti/2- At the time this experi- 
ment was initiated, the '^^Si half-life was believed to be 
in the range of 60 ^ T1/2 ^ 700 yr, and hence a multi- 
year counting experiment was needed to obtain a mea- 
sureable slope. As in other counting experiments, the 
counting rate for ^^Si was continually monitored in the 



same detector against a long-lived comparison standard, 
which in the BNL experiment was ^^Cl (ri/2=301,000 
yr). Since the fractional change in the ^^Cl counting rate 
over the four year duration of the experiment was only 
O(10~^), which was considerably smaller than the overall 
uncertainty of the flnal result, Ti/2(^^Si)=172(4) yr, the 
■^^Cl decay rate was assumed to be constant. Any time 
dependence for ^^Cl beyond the expected statistical fluc- 
tuations was then presumed to arise from various system- 
atic effects, such as drift in the electronics. By computing 
the ratio 32Si/36ci = N{^^Si)/N{^'^C\), these apparatus- 
dependent systematic effects should have largely can- 
celled, and hence this ratio was used to obtain the half-life 
of ^2 Si. On the other hand, barring an accidental can- 
cellation, time-dependent contributions to the '^^Si and 
•^^Cl decay rates themselves would not cancel in the ratio 

32Si/36Cl. 

The BNL data for the ratio ^2 Si/36 ci revealed an un- 
expected annual variation of ^^Si/^^C\ which could not 
be accounted for by the known effects of temperature, 
humidity, or pressure variations in their detector Q. We 
obtained the raw data from the BNL experiment in con- 
junction with an independent effort to apply a new ran- 
domness test 0, [1] to nuclear decays, and the BNL data 
are shown in Fig.[TJ When comparing the results from ex- 
periments on different nuclides, it is convenient to study 
the function U{t) = N(t)/N{0) exp(-t-Ai) rather than 

N{t) itself, since U{t) should be time-independent for all 
nuchdes. For ^^Si, we used A — 4.0299 x 10~^yT~^ from 
Ref. 0. Figure ID exhibits U(t) for the 32si/36ci BNL 
data, along with a plot of 1/i?^, where R is the distance 
between the Earth and the Sun. An annual modulation 
of the 32Si/36ci ratio is clearly evident, as was flrst re- 
ported in Ref. Q. The Pearson correlation coefficient, 
r, between the raw BNL data and l/ii!2 ig 7-=0.52 for 
N=239 data points, which translates to a formal proba- 
bility of 6 X 10~^* that this correlation would arise from 
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two data sets which were uncorrelated. As shown in Fig- 
ure [2l the correlation coeiEcient increases to r=0.65 for 
N—235 data points when a 5 point roUing average is ap- 
pHed. There is also a suggestion in Figs. [1] and of a 
phase shift between and the BNL data, which we 

discuss in greater detail below. 



Hofmallzeil BNL WIM EarHi^un Distance 



J 1.000 




1.03 
1.01 

1.00 I 

% 

CBS " 

o.ee 

0.87 
0.» 



FIG. 1: Plot of U{t) for the raw BNL ^^Si/^GCl ratio along 
with where R is the Earth-Sun distance in units of 

l/(a.u.)^. U{t) is obtained by multiplying each data point 
by exp{+Xt) where A = ln(2)/ri/2 and ri/2 = 172 yr for ^^Si. 
The left axis gives the scale for the normalized U{t), and 
the right axis denotes the values of 1/R^ in units of l/(a.u.)^ 
obtained from the U.S. Naval Observatory (USNO). The frac- 
tional change in ^^Si counting rates between perihelion and 
aphelion is approximately 3 x 10"'^. As noted in the text, 
the correlation coefficient between the BNL data and 1/R^ is 
r=0.52 for A'^=239 points. The formal probability that the in- 
dicated correlation could have arisen from uncorrelated data 
sets is 6 X 10~^*. 



The strong correlation between the BNL decay data 
and the annual modulation of the Earth-Sun distance 
suggests that the ■^^Si/'^^Cl ratio may be responding to 
some influence originating from the Sun. If this is in- 
deed the case, then the effects of this influence would 
be expected to be present in other decays as well. Al- 
though there are hundreds of potentially useful nuclides 
whose half-lives have been measured, the data from many 
of the experiments we examined were generally not use- 
ful, most often because data were not acquired contin- 
uously over sufficiently long time periods. However, we 
were able to obtain the raw data from an experiment 
carried out at the Physikalisch-Technische Bundesand- 
stah (PTB) in Germany 0, [3 

measuring the half-life 
for ^^^Eu, in which ^^^Ra was the long-lived compari- 
son standard. This experiment, which extended over 15 
years, overlapped in time with the BNL experiment for 
approximately 2 years, and exhibited annual fluctuations 
in the ^^^Ra data similar to those seen at BNL. Figure [3] 
exhibits the PTB data as a 5 point rolling average, and 
it is evident from the figure that the PTB data closely 
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FIG. 2: Plot of the 5 point rolling average of U{t) for the 
BNL '^^Si data shown in Figure f . Each data point represents 
the average of 5-points centered on the original datum, which 
serves to smooth short term fluctuations in the "^'^Si/^^Cl ra- 
tio arising from influences other than a possible annual 1/R^ 
variation. As noted in the text, the correlation coefficient be- 
tween the BNL data and I/ii^ is r—0.65 for A'^=235 points. 
The formal probability that the indicated correlation could 
have arisen from uncorrelated data sets is 1 x fO~^^. 



track the annual variation of 1/R^. The Pearson corre- 
lation coefficient r for the data in Fig. [3] is 7-=0.66 for 
A^=1968 data points, corresponding to a formal proba- 
bility of 2 X 10^^^® that this correlation could arise from 
two data sets which were uncorrelated. As in the case of 
the BNL data, there is also a suggestion of a phase shift 
between 1/R^ and the PTB data (see below), although 
this phase shift appears to be smaller than for the BNL 
data. 

Since the BNL and PTB data each exhibit strong cor- 
relations with the annual variation of it is not sur- 
prising that these data correlate with each other. Figure 
|4] exhibits this correlation along with the annual varia- 
tion of 1/R^. The Pearson correlation coefficient for the 
BNL and PTB data is r=0.88 for iV=35 points, which 
corresponds to a formal probability of 4 x 10^^^ that this 
correlation could have arisen from two uncorrelated data 
sets. Moreover, the difference in latitude between BNL 
and PTB, as well as the difference in their climates, ar- 
gues against an explanation of this correlation in terms 
of seasonal variations of climatic conditions such as tem- 
perature, pressure, and humidity etc., which could have 
influenced the respective detection systems. As an ex- 
ample, radon concentrations are known to fluctuate sea- 
sonally, as has been noted in Ref. pL(|], and it was sug- 
gested that the decay of ^^^Rn could lead to a seasonally 
dependent charge distribution on the experimental appa- 
ratus. However, this effect is extremely small given the 
low counting rates that typically arise from radon back- 
ground 



llj, and in any case, the PTB data shown in 



Fig. [3] were corrected for background. 
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FIG. 3: Plot of U{t) for the PTB ^^SRa data along with 
where R is the Earth-Sun distance. See caption to Fig. 1 for 
further details. The fractional change in the ^^^Ra counting 
rates between perihelion and aphelion is approximately 3 x 
10~^. As noted in the text, the correlation coefficient between 
the PTB data and l/R^ is r=0.66 for iV=1968 points. The 
formal probability that the indicated correlation could have 
arisen from uncorrelated data sets is 2 x 10~^*^. Note that 
the la error bars for the PTB data lie within the data points 
themselves. 
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FIG. 4: Correlation between the decay rates of Si at BNL 
and ^^''Ra at PTB. The BNL and PTB data for U{t) have 
been averaged in common weekly bins for purposes of compar- 
ison. The correlation coefficient between the BNL and PTB 
data is r=0.88, which corresponds to a probability of 4 x 10^^^ 
that the BNL/PTB correlation could have arisen from uncor- 
related data sets as a result of statistical fluctuation. Error 
bars are shown for representative BNL data points, and the 
error bars for the PTB data lie within the points themselves. 



The preceding considerations, along with the correla- 
tions evident in Fig.lH suggest that the time-dependence 
of the ^^Si/^^Cl ratio and the ^^^Ra decay rate are being 
modulated by an annually varying flux or field originat- 
ing from the Sun, although they do not specify what 
this flux or fleld might be. The fact that the two decay 
processes are very different (alpha decay for ^^^Ra and 



beta decay for ^^Si) would seem to preclude a common 
mechanism for both. However, recent work by Barrow 
and Shaw 12, 1^ provides an example of a type of the- 
ory in which the Sun could affect both the alpha- and 
beta-decay rates of terrestrial nuclei. In their theory, 
the Sun produces a scalar field (j) which would modulate 
the terrestrial value of the electromagnetic fine structure 
constant aEM- This could, among other effects, lead to 
a seasonal variation in alpha and beta decay rates, both 
of which are sensitive to asM IM- note from Fig. [3] 
that the fractional difference between the ^^^Ra counting 
rates at perihelion and aphelion is w 3 x 10^"^, and this 
would require that the coupling constant ka of to a em 
should be fco, « 3 X 10^. However, this is substantially 
larger than the value ka = (—5.4 ± 5.1) x 10~^ inferred 



from a recent trapped ion experiment |12l . [15| Although 
the specific model of Refs. [l^jE, [l3] would not account 
for the '^^Si and ^^^Ra data quantitatively, variants of 
this model might work. This includes models in which 
separate scalar fields 0i and 02 couple, respectively, to 
ctEM and to the electron-proton mass ratio rrie/mp. 

Another interesting possibility is that terrestrial ra- 
dioactive nuclei are interacting in a novel way with the 
neutrino fiux <i>^ emitted from the interior of the Sun. 
This flux also varies with and the resulting sea- 

sonal modulation of has been observed by Super- 
Kamiokande [l^, 17 1. This possibility is supported by 
the data we report in Ref. [iSj in which we present evi- 
dence for the possible detection of a change in the decay 
rate of ^*Mn during the solar flare of 13 December 2006. 
As noted in Ref. [l8|, the coincidence in time between 
the change in the ^'^Mn counting rate and the solar flare, 
along with other observations, is consistent with a mech- 
anism based on a change in $^ during the solar flare. 

We note that irrespective of the origin of the solar flare 
data, or of the correlations evident in Figs. [l]|4l the ex- 
istence of these effects may explain discrepancies in var- 
ious half-life determinations reported in the literature. 
Examples are ^^Si, ''^''Ti and ^^^Cs, among many others 
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If nuclides such as ^"^Si, ^^Cl, and ^"^^Ra 
respond to changes in the solar neutrino flux due to the 
time-dependence of 1/i?-^, then they can also respond to 
changes in intrinsic solar activity which are known to oc- 
cur over time scales both longer and shorter than one 
year. Thus, depending on when half-life measurements 
were made, and on the specific techniques employed, it is 
possible that some of the half-life discrepancies reported 
in the literature could be reconciled if appropriate data 
on solar activity become available. 

Returning to Fig. 31 we briefly explore the suggestion 
noted above of a possible phase shift of relative 
to both the BNL and PTB data. Although this may 
be an experimental artifact arising from binning effects, 
etc., such a phase shift could also arise from other smaller 
contributions to periodic variations in neutrino flux. Pos- 
sibilities for such contributions were explored in Ref. 17| , 
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where a search was made for short tinie variations in the 
observed flux at Super-Kamiokande arising from either 
the 7.25° incUnation of the solar axis relative to the eclip- 
tic, or from fluctuations in the temperature of the solar 
core. A modulation of the neutrino flux arising from a 
coupling between a neutrino magnetic moment and a lat- 
itudinally inhomogeneous solar magnetic field 2^ could 
also account for a possible phase shift. Although there 
is no compelling evidence at present for such short time 
variations at Super-Kamiokande [Til. \v^ . the statistical 
power of the BNL, PTB, and similar data sets may prove 
to be a useful tool in the search for such effects. Yet an- 
other possible explanation for the apparent phase shift 
could be a seasonally-varying velocity-dependent effect 
similar to that observed by the DAM A/LIBRA collabo- 
ration [Hi. 

In summary, we have presented evidence for a corre- 
lation between changes in nuclear decay rates and the 
Earth-Sun distance. While the mechanism responsible 
for this phenomenon is unknown, theories involving vari- 
ations in fundamental constants could give rise to such 
effects. These results are also consistent with the cor- 
relation between nuclear decay rates and solar activity 
suggested by Jenkins and Fischbach [l^ if the latter ef- 
fect is interpreted as possibly arising from a change in 
the solar neutrino flux. These conclusions can be tested 
in a number of ways. In addition to repeating long-term 
decay measurements on Earth, measurements on radioac- 
tive samples carried aboard spacecraft to other planets 
would be very useful since the sample-Sun distance would 
then vary over a much wider range. The neutrino flux 
hypothesis might also be tested using samples placed in 
the neutrino flux produced by nuclear reactors. 
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